A length-dynamic Tonks gas theory of histone isotherms 
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We find exact solutions to a new one-dimensional (ID) interacting particle theory and apply 
the results to the adsorption and wrapping of polymers (such as DNA) around protein particles 
(such as histones). Each adsorbed protein is represented by a Tonks gas particle. The length 
of each particle is a degree of freedom that represents the degree of DNA wrapping around each 
histone. Thermodynamic quantities are computed as functions of wrapping energy, adsorbed historic 
density, and bulk histone concentration (or chemical potential); their experimental signatures are 
also discussed. Histone density is found to undergo a two-stage adsorption process as a function 
of chemical potential, while the mean coverage by high affinity proteins exhibits a maximum as a 
function of the chemical potential. However, fluctuations in the coverage are concurrently maximal. 
Histone-histone correlation functions are also computed and exhibit rich two length scale behavior. 
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One-dimensional theories in statistical mechanics have 
been successfully applied to numerous biophysical sys- 
tems, including DNA denaturation [1], particle transport 
across biological channels [2, 3], adsorption on ID sub- 
strates [4, 5], and transport [6] along microtubules. In 
this Letter, we study protein-DNA binding and wrap- 
ping by solving a new ID theory of interacting particles 
with dynamically varying particle sizes. DNA- histone 
protein complexes (nucleosomes) play vital roles in com- 
pacting DNA and regulating nucleic acid processing by 
mediating the accessibility by other regulatory proteins 
[7-9]. As shown in Fig. la, DNA can wrap around each 
histone protein complex at most 1.7 times (about 146 
base pairs) [10, 11]. The base pairs that come into con- 
tact with each histone protein defines a footprint which 
we associate with particle lengths. We consider the col- 
lective behavior of the particles mediated only by the 
mutual exclusion of their footprints along a tcnsionlcss 
DNA substrate. We neglect the statistical mechanics 
of the linker DNA between the histone particles as well 
as the non-nearest-neighbor nucleosome-nucleosome in- 
teractions arising from their relative, three-dimensional 
conformations [12]. 

We now derive a generalization of the Tonks- Takahashi 
gas [13, 14] and apply it to DNA-histone structure. Con- 
sider a ID collection of N particles labeled by the po- 
sitions Xi of their left-most edges and confined within 
length L. The first particle is fixed at xq = 0. The mini- 
mum size a of each particle defines the infinite hard core 
repulsive interaction between adjacent particles such that 
\xi — Xi±\ | > a and Na < L. Each particle i also carries 
an additional internal degree of freedom li which corre- 
sponds to its length. The length li may, as we shall see, 
describe the footprint of histone i on the DNA substrate 
(Fig. la). Since thermodynamic properties depend only 
on the interactions among particles, we compute the con- 
figurational partition function 
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FIG. 1: (a) A histone-covered segment of DNA. DNA can par- 
tially wrap nucleosome particles depending on solution ionic 
strength [11, 15]. (b) The ID representation of the footprint 
of the histone- wrapped regions of DNA (thick sections). The 
associated particle lengths can vary from a to a + w. 
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where £* = £*(yi,yi + i) is the maximum possible wrap- 
ping around particle i and e is a mesoscopic free en- 
ergy for unit length extension. This winding-in energy 
consists of the molecular binding enthalpy to the his- 
tone particle, as well as local DNA bending and twist- 
ing energy. If we assume only non-sequence specific 
interactions and uniform extension energies (e = Eq), 
f(y i+ i,yi) = f{y l+ i - yi). The integration limits in (1) 
reflect the particles' impenetrability and renders Z(N, L) 
convolutions of the functions f(yi+i — yi). Upon using 
Laplace transforms in the variable Lq, 
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in the integration limit in (2), we find 
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where /(s) is the Laplace transform of / and 7 e R 
is greater than the real parts of all singularities of the 
integrand. From (3), one readily finds thermodynamic 
quantities such as moments of the particle lengths £i, 
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positions (y™), and particle separations ((yi+i — Di) n )- 

Although these quantities are readily computed using 
specific functions f(s), the resulting sums typically in- 
volve numerical manipulation of extremely large num- 
bers, especially for large N. Thus, it is also useful to de- 
rive from (3), using steepest descents, the leading order 
large N asymptotic approximation Z(N — ► 00, N/L = 
p) ~ e Ar( s *(i/ P -o)+ln /(«*)) ) where s * is the sac idlc point 
defined by largest real root of 
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Thermodynamic limits of, for example, the moments of 
the mean particle lengths (4) become 
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In the histone winding problem (Fig. 1) the unit of 
length will be a single nucleic acid base pair (bp). The 
total length li of each particle corresponds to the arc- 
length of polymer that is wrapped around, and in direct 
contact with, a histone particle. The hard core cut-off 
a depends on details of the three-dimensional arrange- 
ment of adjacent histones, and is roughly (or slightly 
smaller than) the diameter of a histone particle. Only 
for very specific phased orientations of canted histones 
along DNA can the histones be spaced less than about 
a k, 20bp [7, 10-12]. The finite width of the histone lim- 
its "winding-in" lengths £j to either the distance to the 
start of an adjacent particle, Xj+i — Xi — a = yi+\ — Hi, or 
to w ~ 146bp, the maximum winding length correspond- 
ing to 1.7 loops around a histone particle. For example, 
particles zero and one in Fig. la have only one base pair 
of contact (£0 = l\ =0), particles two and three are fully 
wound in (£2 = £3 ~ 146), while particle four is partially 
wound in (0 < £4 < w). Upon imposing the physical 
limits on the particle lengths, 



Xi+\ — Xi — a Xi+i — Xi — a < w 
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With this form of /(s), the integrand in (3) appears to 
have a pole only at s = 0; however, an implicit constraint 
on the number of particles of total length a + £ that can 
fit into L will "induce" poles at s = —Eq. This is seen by 
expanding f N (s) in the integrand of (3), 



Z(N,L) = ^a h * _„,?, ,, v (9) 
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where = (— l) fc (^). For 7 > max{0, — £0} and 
(L — Na — kw) > 0, we close the contour in the left- 
hand s-plane. For L and k such that (L — Na — kw) < 
0, convergence demands that we close the contour in 
the right s-half-plane. Since there are no poles to the 
right of 7, terms with (L — Na — kw) < corre- 
spond to configurations with more particles in L than 
is possible, and do not contribute to the partition func- 
tion. Therefore, we need only sum (9) to k — k* = 
minjint [(L — Na)/w] , TV}. Expanding all terms and ex- 
plicitly evaluating the residues at s — 0, — Eq, we obtain 
the exact expression, 
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where 6 P ee (2(N-l)-p)l/(p\(N-l-p)\). Numerically, 
the asymptotic approximation (6) is accurate to within 
~3% of the exact result (4) provided N > 10. 

Fig. 2a plots the saddle solution s* found from (5). 
In the fixed N ensemble, the probability distribution 
for particle i + 1 to be at position x^i given that 
particle i is at position x\ is readily computed in the 
asymptotic limit [L ^> A ^> a), g^ 2 \xi + \ — Xi\N) ~ 

Z(N,L~x)/ J a A Z(N,L-x)dx = s*e~ s '( x - a \ This result 
implies that the adjacent particle is statistically confined 
to within x < 1/s* [16]. At low number densities, adja- 
cent histones are spaced far apart and s* ~ p/(l — po)- 
For attractive lengthening interactions (eq C 0), a sharp 
increase in s* occurs near p ~ l/(w + a) signaling a 
partial confinement of hard rod particles of roughly size 
w + a. At extremely high densities, p ~ a" 1 , particles 
are compressed at the expense of unwinding, and s* in- 
creases further as s* ~ 2p/ (1 — pa) — £q/2. Fig. 2b shows 
the mean winding-in length (normalized by w), found 
from (6). At low densities and strong attractive binding, 
the maximal winding in length w' 1 (£) ~ 1 is approached, 
while at high densities, the winding-in length is restricted 
by nearest neighbors and (£) « 1/p. 
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FIG. 2: Thermodynamic properties for a = 20, u> = 146, 
and various eo. (a) The root s* of (5) that defines Z(N — > 
oo, N/L = p) and determines the next-neighbor correlation 



function '(xi+i — Xi\N) = s*e~ 



3 (x — a) 



(b) Winding-in 



lengths as a function reduced density pw. (c) Num- 

ber density as a function of chemical potential p for various 
binding affinities eqw. An intermediate density plateau at 
{p)a ~ a/(w + a) arises for high affinity histones due to close 
packing of fully wound-in particles of length w + a. (d) Mean 
coverage as a function of p. Curves correspond to the leg- 
end used in (b). Maxima arise for wound-in, high affinity 
particles spaced an average distance w + a apart. The high 
compression (p — > oo) limit forces unwinding but packs his- 
tones at intervals of the linker length a for a coverage fraction 
(9) a l/a. 



Histone-DNA affinity and competition experiments 
however, are performed by exposing DNA to fixed bulk 
histone concentrations [11, 17]. When the mean bound 
histonc number is determined by the bulk histone chem- 
ical potential (thus not necessarily large) , we employ the 
grand partition function found using the exact expression 
(3) or (10) in 3(A, L) = J2n"=i X * Z ( N > l )- Thc fugacity 



A = 



60 takes into account the bulk histonc chemical 



potential, and the binding energy £o of a single base pair. 
N* = int{i/a} is the maximum number of particles that 
can fit into length L. In thc L/a = oo, N* = oo limit, 
we find 
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which can be evaluated using (8) to yield 3 (A, L — ► oo) = 
(p)e s+L , where the mean density (p) = L _1 cL ln3 is ex- 
plicitly 



(P) = 



S+(S+ + £p) 

2s+ + e + Ae~ s + a (a - (w + a)e-0++ E °)<") 



, (12) 



s+ being the largest real root of 1 — A/(s + )e~ s+a = 0. 
The mean fraction of DNA base pairs covered by contacts 
with histones is then found from 
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The mean density and coverage are plotted in Figs. 



2c, d. For fx 



3, (p) ~ s+ ~ e- £ »+"(l 



0/eo 



since this limit corresponds to infinitely dilute bulk hi- 
stone concentration. Densities of adsorbed particles in- 
crease with bulk histonc chemical potential. For higher 
affinity histones, these increases occur earlier and plateau 
at a value corresponding to close packing of fully wound- 
in particles, or (p)a ~ a/(w + a) ks 0.12. As p and the 
density further increase, particles compress and unwind 
each other until the density approaches maximal pack- 
ing at (p(p — ► oo)) ~ \/a - 2/ '(pa) + 0(p~ 2 ln/i). This 
crossover behavior is seen only for high affinity (eq <C 0) 
histones. Thc transition from monotonic to the two-stage 
density behavior may be observable by tuning thc non- 
specific histone binding energies eo by e.g. varying ionic 
strength [11, 15]. 

DNA accessibility (and hence occlusion by histone par- 
ticles) and positioning is thought to be a major determi- 
nant of nucleic acid processing by other regulatory pro- 
teins [7-9, 11]. The total coverage (8) defined as the mean 
fraction of base pairs in contact with histones, is shown in 
Fig. 2d. For p — > — oo, there are few particles present to 
adsorb onto the DNA substrate and (0) ~ e ~ e ° +fl -» 0. 
An increase in mean coverage results from increasing p 
and the number of bound histones. The theoretical max- 
imum coverage w/(w + a) m 0.88 (thin dashed line) cor- 
responds to close packing of fully wound-in particles and 
is approached only for high affinity histones. The mini- 
mum uncovered fraction a/ (w + a) results from the linker 
DNA of minimum length a ~ 20 joining two adjacent mi- 
crosomes. If pt is further increased, the particles squeeze 
on each other until they unwrap to the point where only 
a single base pair remains in contact for each histonc. 
In this limit, the mean coverage (9(p — > oo)) ~ l/a + 
(1 — 2/a)/p + 0(p~ 2 \np) ss 0.05 while the histones are 
spaced at their maximal densities (p) w l/a. The vari- 
ance in coverage, var(0) = L~ 1 d(L~ 1 dE/deo + S)/d£o, 
gives a standard deviation in coverage proportional to the 
mean coverage, \J var(0) = y/2(0). Despite the high cov- 
erage at intermediate chemical potentials, fluctuations in 
this regime are also maximal, suggesting a dynamically 
controlled DNA accessibility mediated by histone-histone 
interactions. 

In the grand ensemble, the correlation function analo- 
gous to g^^x^-x^N) is (p)g^(x\p) = 3(A, x)S(X,L— 
x)/S(A, L) which describes the probability distribution 
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that given a particle at the origin, any other histone ex- 
ists at x. When x ^> 1 / (p) and N* < oo, the exact 
expression (10) must be used to compute H(A, x). Fig. 3 
shows g^ (x\p) computed for various values of £o, fJ>- Fi- 
nite particle size anticorrelations give rise to oscillations 
at two length scales (see Fig. 3 caption) depending on 
affinity eo and density (p(p)). 




FIG. 3: The correlation g^ 2 \x\p) exhibits properties of both 
length scales a and w depending on eo and p. For high affinity 
and low densities {e.g. eow — — 30, p = —5), the density 
distribution is similar to that of a hard rod Tonks gas of length 
a + w . Upon increasing the density (eqw — — 30, p = —2), 
features associated with both length scales arise as partial 
unwinding occurs, exposing the hard core repulsion of size a. 
Conditions under which particles weakly wind in (sow — —10) 
exhibit behavior attributed to hard rods of length a. However, 
even at lower affinities (eo^ = — 8), behavior approximating 
that of rods of length w + a can be recovered if densities are 
made sufficiently low (p = —15). 

The proposed model considers only the adjacent his- 
tone exclusion interactions mediated by nonoverlapping 
footprints and neglects sequence specific and nucleosome- 
nucleosome interactions arising in compact, 3D chro- 
matin structure. Nonetheless, our theory can be solved 
exactly with uniform wrapping energies eo to give rea- 
sonable results for winding-in lengths (£) and histone- 
histone correlations g( 2 \x\N) and g^(x\p). The model 
predicts a two-state adsorption process and a maximum 
in the DNA coverage fraction (9) as a function of bulk 
histone concentration and binding affinity. However, fluc- 



tuations in the coverage are also concurrently maximal 
during the peak in mean coverage, indicating that ther- 
mal effects can nonetheless provide dynamic access to 
highly covered DNA. These predictions can be tested ex- 
perimentally by varying bulk histone coverage and ionic 
strength (although this would also affect DNA bendabil- 
ity and eo)[ll, 15], provided short enough DNA segments 
are used such that large scale 3D structures do not arise 
and non- nearest-neighbor interactions remain irrelevant. 

The model can be readily generalized to include the 
effects of linker DNA twist, externally applied tension 
[18], and relative histone orientation [12] on the cut-off a. 
Specific sequences, and their effects on local bendability, 
twistability, and histone affinity has also been found to be 
important in vitro [11] and can be treated within a similar 
framework, although completely generalizing our model 
to include specific sequences [11, 17] (spatial dependence 
of e(x)), would require computational approaches. 

Finally, we can extend the Tonks- Takahashi model by 
removing the maximal winding length constraint (w — > 
oo). The resulting model can then be applied to problems 
of ID nucleation and polymerization. Particles would 
correspond to a nucleation domains that cannot coa- 
lesce due to perhaps growth asymmetry (as in actin fil- 
aments) or incommensurability with the underlying sub- 
strate [5, 19]. The "winding-in" or polymerization length 
£ is now the degree of polymerization in that domain. 
Since £i are limited only by the position of adjacent do- 
mains, I* = x%+\ — Xi — a and f(s) = [s(s + eo)] -1 - Re- 
sults analogous to those obtained for the histone problem 
are found by setting w = oo in (12-13). The qualita- 
tive behavior of £(eq) can be found in closed form in the 
N = oo,a = limit: p n {£%) ~ 2™n!/(z + 2 + Vz 2 + 4)", 
where z = eq(1 — pa)/ p. This result also implies 
{£ 2 ) — (£i) 2 = {£i) 2 . The behaviors for w = oo are qual- 
itatively different from those shown in Fig. 2. Namely, 
the adsorption isotherm is monotonic and the coverage 
has a minimum at intermediate p when eo < and van- 
ishes for p — > — oo if eo > 0. Such differences highlight 
the importance of the additional length scale w unique 
to histone binding. 
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